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Abstract. In a standard cold dark matter (CDM) cosmology, microhalos at the CDM cutoff 
scale are the first and smallest objects expected to form in the universe. Here we present 
I results of high resolution simulations of three representative roughly Earth- mass microhalos in 

2 order to determine their inner density profile. We find that CDM microhalos in simulations 

without a cutoff in the power spectrum roughly follow the NFW density profile, just like 
the much larger CDM halos on galaxy and galaxy cluster scales. But having a cutoff in 
the initial power spectrum at a typical neutralino free streaming scale of 1O~^M0 makes 

CN their inner density profiles considerably steeper, i.e. p oc 7--(i-3-i-4)^ good agreement 

with the results from Ishiyama et al. (2010). An extrapolation of the halo and subhalo 
mass functions down to the cutoff scale indicates that microhalos are extremely abundant 
throughout the present day dark matter distribution and might contribute significantly to 
indirect dark matter detection signals. Assuming a transition from a NFW to a steeper inner 
profile (p oc r~^'^) two orders of magnitude above the cutoff scale, the total boost factor for a 
Milky Way sized dark matter halo increases from about 3.5 to 4. We further find that CDM 
microhalo concentrations are consistent with the Bullock et al. (2001) model and clearly 
^ rule out simplistic power law models for the mass dependence of concentrations and subhalo 

^ annihilation, which would erroneously lead to very large boost factors (a few hundred for 

galaxy halos and over 1000 for clusters). 
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1 Introduction 

The standard model of cosmology is characterised by a hierarchical bottom-up formation of 
structure [1], in which the size of the smallest objects is set by the free streaming length of 
the dark matter particle. For a standard supersymmetric cold dark matter (CDM) candidate 
with a mass of 100 - 1000 GeV, free streaming leads to a cutoff in the matter power spectrum 
at a scale in a range of about 10~^^ to 10~^ Mq [e.g. 2-9]. Such microhalos, and especially 
their inner regions, might be dense enough to survive tidal disruption up to present time and 
they are therefore potential contributors to indirect dark matter detection signals [10-19]. 

The density profiles of CDM halos from dwarf galaxy to galaxy cluster scales have been 
studied extensively with many high resolution cosmological simulations [e.g. 20-24] . All CDM 
halos on theses scales roughly follow the universal NFW profile [25] . The NFW function has 
only one free parameter, the scale radius or concentration c = J'viriai/'^s, which is directly 
related to the initial mass fluctuation spectrum and the typical halo formation times [26-29] . 
Unlike these well studied massive CDM halos, the first, smallest microhalos do not form 
through hierarchical merging and they do not contain any subhalos. Their internal structure 
might be different from that of cluster- or galaxy-sized dark matter halos. The first numerical 
simulations of microhalo formation [14] did not have the resolution to resolve the central 
regions of microhalos, and in the well resolved radial range the microhalos showed similar 
profiles as those found in large CDM halos. Recently, Ishiyama et al. simulated microhalo 
formation with 20 times better mass resolution and surprisingly steep inner slopes were found 
for halos near the cutoff scale [30]. Simulated without the cutoff in the initial power spectrum, 
the same microhalos turned out to have shallower NFW-like inner profiles. In this work we 
simulate the formation of microhalos with a mass resolutions as low as nip = 1.46 x 10^^'^ 
Mq, 64 times higher than Ishiyama et al. [30], and we are able to largely reproduce their 
surprising results about the steeper density profiles near the cutoff scale. Substructure in 
the form of microhalos contributes significantly to the dark matter annihilation signals from 
all larger scale dark matter objects (halos, subhalos, galactic center). Here we show that the 
steeper mircohalo density profiles might increase the substructure boost factor of a galaxy 
sized CDM halo by about 5 to 12 per cent. 
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2 A^-body Simulations 



Numerical simulations have been performed with the parallel treecode pkdgrav, written 
by Joachim Stadel and Thomas Quinn [31]. Initial conditions are generated with a parallel 
version of the GRAFIC package [32] and are based on the cosmological parameters taken 
from WMAP7 [33]: ag = 0.8, h = 0.7, = 0.227, = 0.046, JIa = 0.727 and = 0.961. 
Our simulations start at an initial redshift of z^c = 500 and run until Z{ = 31, have a force 
softening length e of 1/50 times the mean inter-particle separation and a fixed time-stepping 
with A^steps = lO'OOO. 

To begin with, we have run two cosmological boxes of 30 pc and a total particle number 
of N = 512^, both having identical random seeds. One of the simulations contains an 
exponential cutoff in the matter power spectrum corresponding to a CDM particle mass of 
~ 1 TeV [5]. We use the friends-of- friends algorithm [34, 35] with a linking lenght of b = 0.2 
to identify the halos at Zf = 31. Then we select the three largest halos and rerun them at 8 
times better mass resolution, again with and without cutoff. These refined simulations have 
a mass resolution of nip = 9.3 x 10^^'^ M0 and a softening length of e ~ 5 x 10^^ pc, which 
is the same as in the Ishiyama et al. (2010) simulations [30] . 

In order to test the very inner structure of microhalos right at the cutoff scale, we 
rerun one of the three halos once more with further increased resolution, this time only 
for the case with the cutoff in the power spectrum. This ultrahigh resolution microhalo 
contains approximately 33 x 10^ particles within r2oo with a simulated particle mass of 
nip = 1.46 X 10""^^ Mq (64 times better mass resolution than in [30]) and a softening length 
of e ~ 1.25 X 10^^ pc. For what follows, we will name the two levels of resolution Level 
1 and Level 2 respectively. To check for numerical convergence we did rerun our Level 1 
simulation with different starting redshift (zic = I'OOO), force resolution (e/2) and time- 
stepping (A'steps = 5'000), see Section 3 and the last panel in Figure 2. 

Figure 1 shows the projected density (top-row) and phase-space (bottom-row) density 
maps of the Level 2 halo at four different redshifts {z = 31, 33, 39, 43). At z = 31, there are 
almost no subhalos, but there is abundant substructures in the form of caustics, similar to 
those in the classical self-similar secondary infall model [36, 37]. Such caustics could increase 
the dark matter annihilation rate by contributing to the total annihilation luminosity (see 
Section 4) [38-40]. To estimate this contribution for the present Level 2 halo, we compare the 
signal coming from each individual simulation particle within the halo (Li ~ Pinii), with 
the one coming from when binning up the particles in spherical shells and integrating over 
the volume (L2 ~ / p(?')bm'^^)- ^^"^ ratio L1/L2 to be 3.5 at z = 43, and decreasing 
to roughly 2.3 at 2: = 31. By smoothing out these caustics (i.e. increasing the number of 
considered neighbours from 32 to I'OOO when computing the local densities), L1/L2 shrinks to 
1.5 at z = 31. This suggest that at 2; = 31 the microhalo annihilation luminosity is increased 
by a factor of about 1.5 ~ 2.3/1.5 due to caustics and by another factor of 1.5 because of other 
departures from spherical symmetry. It is possible that the actual contribution from caustics 
is significantly larger than that, as we only start to resolve caustics in our highest resolution 
simulations and more detailed future simulations might resolve more and sharper caustics. 
However, we do see that relative contribution of caustics to the microhalo annihilation signal 
decreases with time and it is difficult to estimate how large it will be at 2; = 0. Therefore we 
decided to ignore their potential contribution when calculating the boost factors at z=0 in 
Section 4. 
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Figure 1. Top: Projected density maps of the high resolution microhalo, at redshifts (from left to 
right): z — 31, z = 33, z — 39, z — 43. Bottom: corresponding phase space density maps, calculated 
with Enbid. In all panels, (phase space) densities are represented in a logarithmic color map: bright 
regions refer to a high (phase space) density whereas dark regions refer to a low (phase space) density. 



3 Microhalo density profiles 

In order to quantify the detectability of microhalos (e.g. [30]) and the survival probability 
as galactic subhalos [10-16], one needs to know the form of their density profiles and their 
typical concentrations (we use the concentration definition C200 = '''200/rs, where r2oo is the 
radius enclosing 200 times the critical density of the Universe). 

3.1 The inner density profiles 

Figure 2 shows the spherically averaged density profiles of the three largest Level 1 microhalos 
(see Table 1 for details) with (red triangles) and without (black squares) cutoff. It is evident 
that the inner slope of the models with cutoff is substantially steeper than without the cutoff. 
In order to quantify the central structure of the three halos in the cutoff model, we have fitted 
the density profiles with the following conventional parametrisation [e.g. 25]: 

p{r) = ^ (3.1) 

Using the Levenberg & Marquardt method, we obtained the following values for the slope of 
the inner profile via a minimisation: 

"[Halo 1] - 1-4, a[Halo2] ~ 1-3, a[Halo 3] - 1 •4- (3.2) 

The fits are shown by the red solid lines in Fig. 2 (panels 1-3). On the other hand, a NFW 
profile (a = 1) provides a reasonable fit for the microhalo densities in the model without a 
cutoff (black solid lines). 
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Figure 2. Panels 1-3: Spherically averaged density profiles of the three largest collapsed microhalos 
at z = 31, with (red triangles) and without (black squares) cutoff. The red solid lines refer to the 
best fit according to Eq. (3.1) with a = 1.4 (Halo 1 & Halo 3) and a — 1.3 (Halo 2), the black solid 
lines refer to a NFW fit respectively. The radial distance is plotted in physical units, densities in 
units of pcrit at z = 0. Panel 4: Density residuals between the Level 1 run and three convergence test 
simulations, each varying one simulation parameter. 



In order to exclude a systematic computational error leading to such steep profiles, we 
have rerun the Level 1 simulation in three variations to verify numerical convergence: (i) 
with half as much time-steps (A'^gtops = 5000), (ii) with a starting redshift of z\c = I'OOO 
and (iii) with a reduced gravitational softening length (e/2). The resulting density profiles 
are represented as residuals {pi/ p — 1) in the last panel of Figure 2. In all three halos, the 
deviations are negligible for the A^^^teps = S'GOO (blue solid lines) and the e/2 (green solid 
lines) simulation. The run with increased starting redshift (red solid lines) shows the largest 
fluctuations, most probably due to the fact that force errors become significant in the nearly 
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Table 1. Halo parameters of the Level 1 simulation at redshift z — il. M200 and r2oo are measured 
as 200 times the critical density, a is the inner density slope of the measured density profile (see Eq. 
(3.1)), a = 1 corresponds to the NFW profile. Distances are given in physical units. 
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uniform, very high redshift matter density field. However, even for the latter case, differences 
are only of the order of 20 per cent. We conclude that our results are numerically robust, and 
that the origin of such steep density cusps must indeed be physical. The different inner profiles 
of microhalos near the cutoff scale is most likely related to their different, non- hierarchical 
formation, but we are currently unable to explain this surprising result. 

Halo 2 was resimulated at the Level 2 refinement, i.e. with 64 times better mass resolu- 
tion. This allows us to resolve the very inner part of the density profile, down to a resolved 
radius of ~ 5 x 10"^ pc. Even at these very small scales we do not observe any flattening or 
deviation from a a = 1.3 cusp, as indicated by the dashed black line in Figure 3. Further 
the evolution of the density profile from z = 43 to 2: = 31 is shown. The inner slope of the 
microhalo remains constant after about z = 39. 

3.2 Microhalo concentrations 

Dark matter halo concentrations are directly related to the halo formation times, which can 
be calculated from the mass fluctuation spectrum a{M) (see e.g. [26-29]). On microhalo 
scales, (t(M) decreases only very slowly with mass, which leads to nearly constant typical 
concentrations. The microhalo simulations in [14] used a box size of 3 comoving kpc which 
contained a small, very high resolution region where microhalos (M ~ 10^^ Mq) were re- 
solved. The surrounding region with lower resolution contained some larger halos (M ~ Mq) 
which indeed showed similar concentrations as the microhalos. The nearly constant a{M) 
also means that the scale of the small simulation boxes (or refinement regions) used to resolve 
microhalos enter the non- linear regime soon after the typical microhalo formation time. This 
is the reason why microhalo simulations usually end around z = 30 and usually resolve only 
a small sample of objects. Finite box size simulations also exclude large scale fluctuations, 
which are significant in this mass and redshift regime. A 30 pc box for example, as used in 
this work and in Ishiyama et al. (2010) [30], could lower the typical formation redshift by 
about a factor of two (see Figure 3 in [30] ) , which would lower the microhalo concentrations 
by the same factor. On the other hand, the microhalos which exist in our final snapshot 
(z = 31) represent an earlier forming, more concentrated subset of the full microhalo sample. 

Besides the limitations mentioned above, we will still use our sample of microhalo con- 
centrations for a rough comparison with CDM mass - concentration relations. At z = 31 and 
M = lO"'^ Mq, the Bullock et. al (2001) model ([26], in the normalisation of [27]) predicts a 
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mean concentration cnfw,200 = 1-85 with a halo-to-halo one sigma scatter of A(logc) ~ 0.18 
at a given mass. A more recent, refined version of the Bullock et al. (2001) model [28] 
predicts cnfw,200 = 2.43 and the same scatter. Both models are consitent with the cnfw,200 
values found here (see Table 1) and also with the microhalo samples from earlier microhalo 
simulations [14, 30]. 

In the simulations with a cutoff the profiles do not follow NFW and we used steeper 
functions {a = 1.3, 1.4). This shifts the best fit scale radius further out and therefore 
reduces the concentration parameter compared to the NFW case. Keeping the radius ry^^^ 
of the maximum circular velocity fixed, one can convert the concentrations between an NFW 
profile and variations of it [e.g. 59, 60]. For the case of a = 1.3 and a = 1.4 respectively, the 
conversion is given by 

Cl.3 = CNFw/1-42, Ci.4 = CNFw/1-67. (3.3) 

The cnfw values obtained from this conversion (see Table 1) are higher than those of the 
halos in the no-cutoff model. Apparently the cutoff increases both the microhalo central 
densities and their concentrations. 

After halo formation the typical scale radii of halos of a given mass are roughly constant, 
while r2oo and C200 increase proportional to the expansion factor a = 1/(1 -|- z) [e.g. 26]. 
Assuming this average scaling relation applies also to our individual microhalos, the z=0 
concentration estimates for Halo 1-3 in the cutoff model is C200 = 97.0, 93.1 and 67.2. These 
are again in line with mean concentrations predicted by the models at z=0: cnfw,200 = 59.2 
(Bullock et al. (2001)) and cnfw,200 = 77.8 (Maccio et al. (2008)). 

A simple power law approximation to the Bullock et al. model is often used on galactic 
halo scales, cnfw,200 — 8.45 x (M/IO^^Mq)"*^-^^. Extrapolating this simple approximation 
all the way down to microhalo scales (10"'' M©) give typical microhalo concentrations of 
around 1040 at z=0. Even with our small sample of z=31 microhalos (or with any of the 
earlier microhalo samples from [14, 30]) one can clearly rule out such a simple power law 
concentration-mass relation. 

4 Small scale structure contributions to the DM annihilation signal 

Dark matter annihilation rates are proportional to the density squared, so any small scale 
clumpines contributes to indirect dark matter detection signals. The total annihilation lumi- 
nosity L from a dark matter halo consists of a main halo signal L from the smooth spherically 
averaged dark matter distribution and from a substructure contribution Lgubi often expressed 
as a halo boost factor B{M): 

L{M) = L{M) + Lsub(M) = L(M) + B{M)L[M) = L(M)(1 + B{M)), 

Substructure boost factors play a significant role in the prediction of indirect detection 
signal, for example when deciding if nearby small dwarf galaxies or local galaxy clusters are 
the most promising targets [41]. Large A^-body simulations have shown that the numerically 
resolved subhalos in galactic halo have L^nh — LsnhiM), i.e. B{M) ~ 1, and that significant 
contributions to B(M) from unresolved smaller subhalos gives boost factors of a few up 
to about ten, depending on details of the extrapolations to microhalo scales [24, 42]. On 
the other hand, simplistic models based on power laws for c(M) or Lsub(> M) give much 
larger boost factors of a few hundred for galactic halos [43], see [44] for a recent review and 
discussion. The different boost values are mostly caused by very different assumptions about 
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Figure 3. Redshift evolution of the spherically averaged density profile in the Level 2 refinement of 
Halo 2. The black dashed line indicates the inner slope of a = 1.3. Distance and density are plotted 
in physical units. 



microhalo concentrations. In the following we will show that only the lower estimates (e.g. 
for galactic halos B{IQ^'^Mq) < 20) are consistent with simulated microhalo concentrations. 

The boost factor of a dark matter halo of mass M is given by the integrated subhalo 
contribution (form the CDM free streaming cutoff mass mo up to the main halo mass M) 
divided by the luminosity of the main halo L [e.g. 44-49]: 



where dN/ dm is the subhalo mass function and a its slope. The normalisation of the subhalo 
mass function, A, is set by requiring a fraction of the main halo mass to be in subhalos. For 
a = —1.9 it is obtained to be A ~ 0.03 [e.g. 50, 51]. The total annihilation luminosity coming 
from a single halo of mass M is given by 



where c is the concentration parameter and /(c) is fixed by the enclosed mass (and is therefore 
dependent on the exact form of the density profile) [e.g. 52]. For a two parameter density 




(4.1) 



L(M) = AttM 
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(4.2) 
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profile as in Eq. (3.1), the general form of /(c) is given by 



/(c) = 1, X 2Fi(3 - a; 3 - a; 4 - a; -c), (4.3) 

\6 — a/ 

where 2F1 is the hypergeometric function. In the case of a NFW profile (q = 1, 7 = 3), Eq. 
(4.3) reduces to /(c) = log(l + c) — c/(l + c). To evaluate B{M) a one needs to know the 
form of the subhalo density profiles and also thier concentrations. 

4.1 Moderate microhalo concentrations rule out very large boosts from power- 
law extrapolations 

For the subhalo mass function the extrapolation of a simple power law is consistent with 
theoretical models (e.g. [53] and with microhalo simulations [14, 54]). The same is not 
the case for the concentration-mass relation c(M): models and microhalo simulations are 
consistent with concentrations related to halo formation times and therefore to the mass 
fluctuation spectrum a{M) and to the power spectrum P{k), which are not power laws (see 
section 3.2). In realistic models one combines a power law subhalo abundance with a non- 
power law concentration model and gets a total subhalo luminosity Lsub("^o)) which is not a 
simple power law of the minimum subhalo mass mo- The very large boost factors (over 200 
for galaxy halos [43] and over 700 for clusters [55]) obtained from power law extrapolations 
are ruled out by these and earlier microhalo simulations [14, 30] and by the typical halo 
formation times expected from the CDM mass fluctuation spectrum a{M). 

We use the concentration - mass relation of Maccio et al. (2008) [28], which is consistent 
with microhalo simulation results (see section 3.2). Assuming universal NFW profiles for all 
subhalos down to mo = 10""^ M© (and up to M = 10^^ M0) gives the moderate value of 
B{M) = 3.56, as indicated by the black solid line in Figure 4 (which is in agreement with 
the results of other authors [e.g. 14, 45, 46]). This will be our reference value derived from 
a simple, plausible subhalo model. Slightly different models are possible given our current 
understanding of structure formation and could include the following variations which change 
the boosts by small factors: 

• Instead of NFW one could use the Einasto profile (with a fixed ue = 0.17) instead. 
It fits density profiles of large CDM halos better than NFW [e.g. 22, 56], it is denser 
than NFW in the inner halo (around 0.003 to 0.1 r^) and normalised to the same outer 
profile and concentration gives about 1.4 times larger annihilation fluxes[57, 58]. 

• We neglect halo-to- halo scatter in the concentration parameter [26, 28] and the profile 
shape [22, 23] (i.e. in the second parameter a^;, when the Einasto form is used), which 
could lead to a small increase in the boost factors. 

• The subhalo mass functions might be slightly steeper, even the critical value of a = —2.0 
(i.e. constant mass in subhalos per mass decade) is consistent with some simulated 
CDM halos [42, 55]. Changing a from -1.9 to -2.0 in our reference NFW based model 
would increase the boost by about a factor of 6 {B{M) = 21.7). 

• We neglect caustics in microhalos. Their contribution to the z = annihilation signal 
is uncertain: on the one hand the effect our simulations show at 2 = 31 might increase 
with numerical resolution and on the other hand the evolution from z = 31 until present 
time is unclear. However, since their moderate contribution at z = 31 appears to decay 
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Figure 4. The annihilation himinosity boost factor extrapolated down to a cutoff mass of mo = 
1O^^M0. The black solid line shows the boost factor results under the assumption of NFW density 
profile {a — 1) in combination with a concentration-mass relation as suggested by [28]. The dashed 
colored lines show the additional contribution to the boost factor when the inner density profile 
of the smallest microhalos has slope of a = 1.3 (red), 1.4 (blue) and transition scales of lO^^M© 
(short-dashed) and lO^^M© (long-dashed) respectively. 



quickly with time (see Figure 1 and Section 2), we do not expect them to affect z = 
boost factors significantly. 

• Steeper inner density profiles for subhalos near the cutoff scale, as found in Ishiyama 
et al. (2010) [30] and in the simulations presented here, increase the total annihilation 
boost factor moderately, as described and quantified in the next section. 

4.2 Steeper microhalo inner density profiles increase boost factors moderately 

Microhalos near the cutoff scale have steep central density cusps with a slope of a = 1.3 — 1.4. 
In order to quantify the contribution of such objects to the total annihilation boost, we solve 
Eq. (4.1) numerically with the boundary condition B{mQ) = 0. In addition, we introduce 
a transition mass Mj = {10~^, 10~^} Mq as a sharp transition scale: above Mt the density 
profiles are assumed to be NFW like, below Mt they are assumed to have the functional 
form found in Section 3. Fixing the cutoff scale to be mo = 10^^ M0, the results are 
shown in Figure 4. The dashed colored lines show the increase of B{Ai) when the smallest 
objects possess a slope of the inner density profile of a = 1.3 (red) and a = 1.4 (blue), with 
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transition scales of Mt = 1O~^M0 (short-dashed) and Mt = 10~^ Mq (long-dashed). We 
have incorporated the concentration correction (see Eq. (3.3)) in the calculations of B{M). 
It is evident from Figure 4 that such steep microhalo cusps increase the total annihiliation 
boost for a Milky Way sized halo between 5 to 12 per cent, strongly dependent on the inner 
slope, the transition mass and the cutoff scale. 

In the analysis above, we have assumed that the slope of the inner density profile stays 
unperturbed between z = 2>1 and z = 0. Unfortunately, current state of the art simulations 
do not allow to evolve a microhalo in a realistic environment to the present epoch. However, 
there is good reason to assume that the inner density profile of microhalos stays unchanged 
until z = 0: on the one hand, it was shown in earlier works, that the inner density profile of 
larger virialised halos is practically unaffected by mass accretion [e.g. 61, 62]. Furthermore 
recent results show that in mergers of collisionless dark matter halos the cusp of the steeper 
progenitor is preserved [e.g. 63, 64]. 

5 Summary & Conclusion 

We have used ultra high resolution A^-body simulations to examine the slope of the inner 
density profile of microhalos and the associated implications for the total annihilation lumi- 
nosity boost factor of a galaxy sized dark matter halo. Our results can be summarised as 
follows: 

• When introducing an exponential cutoff in the matter power spectrum at a scale that 
corresponds to the size of a typical microhalo, we found that the density profile of these 
objects is accurately described by a NFW-like fitting function (see Eq. (3.1)) with an 
inner slope of q = 1.3 — 1.4. Without the cutoff, the inner structure is found to be very 
similar to the one of larger, i.e. galaxy- and cluster-sized, dark matter halos, and are 
reasonably well fitted by the usual NEW form. 

• Assuming that one (or two) order of magnitude above the cutoff scale the inner slope 
of the density profile rapidly changes from a = 1 to a = 1.4, we find that the total 
annihilation boost factor of a Milky Way sized halo is increased by 5.8 (11.9) per cent. 
If the inner slope is q = 1.3, the increase is found to be 4.7 (9.6) per cent. 

• The measured concentrations of microhalos agree with the concentration - mass relation 
of Maccio et al. (2008) [28] (and therefore also with the Bullock et al. (2001) model 
[26]). Power law models for the concentration - mass relation, which would boost the 
total annihilation luminosity by two to three orders of magnitude, are therefore clearly 
ruled out. 

Larger microhalos simulations, like the recently presented trillion particle simulation [65], 
will be needed to reduce some remaining uncertainties, especially about the transition mass 
scale, at which the inner slope of the density profile significantly increases, the scatter in 
microhalo profiles and the exact mean microhalo concentrations and scatter. More work is 
also required to understand the relations between inner profiles, halo formation histories and 
cutoff scales. 

Acknowledgments 

We would like to thank Mike Kuhlen, Miguel Sanchez-Conde and Neal Dalai for helpful dis- 
cussions and comments in the final stages of this work. Numerical simulations were performed 



- 10 - 



on the SuperMUC at the LRZ in Garching, on the zBox3 at the Institute for Theoretical 
Physics in Zurich and on the Rosa cluster at the CSCS in Lugano.This work was supported 
by the SNF. 



References 

[1] P. J. E. Peebles, Large-scale background temperature and mass fluctuations due to 
scale-invariant primeval perturbations, 1982, ApJ, 263, LI 

[2] A. Loeb & M. Zaldarriaga, Small-scale power spectrum of cold dark matter, 2005, Phys. Rev. 
D, 71, 103520 

[3] E. Bertschinger, Effects of cold dark matter decoupling and pair annihilation on cosmological 
perturbations, 2006, Phys. Rev. D, 74, 063509 

[4] G. Bertone., D. Hooper, J. Silk, Particle dark matter: evidence, candidates and constraints, 
2005, Physics Reports, 405, 279 

[5] A. M. Green, S. Hofmann, D. J. Schwarz, The power spectrum of SUSY-CDM on subgalactic 
scales, 2004, MNRAS, 353, L23 

[6] S. Profumo, K. Sigurdson, M. Kamionkowski, What Mass Are the Smallest Protohalosf, 2006, 
Phys. Rev. Lett., 97, 031301 

[7] T. Bringmann, Particle models and the small-scale structure of dark matter, 2009, New Journal 

of Physics, 11, 105027 

[8] S. Hofmann, D. J. Schwarz, H. Stockcr, Damping scales of neutralino cold dark matter, 2001, 
Phys. Rev. D, 64, 083507 

[9] R. Kolb, L I. Tkachev, Femtolensing and Picolensing by Axion Miniclusters, 1996, ApJ, 460, 
L25 

[10] H. Zhao, D. Hooper, G. W. Angus, J. E. Taylor & J. Silk, Tidal Disruption of the First Dark 
Microhalos, 2007, ApJ, 654, 697 

[11] T. Goerdt, O. Y. Gnedin, B. Moore, J. Diemand & J. Stadel, The survival and disruption of 
cold dark matter microhaloes: implications for direct and indirect detection experiments, 2007, 
MNRAS, 375, 191 

[12] V. Berezinsky, V. Dokuchaev, Y. Eroshenko, Small-scale clumps in the galactic halo and dark 
matter annihilation, 2003, Phys. Rev. D, 68, 103003 

[13] V. Berezinsky, V. Dokuchaev, Y. Eroshenko, Small-scale dark matter clumps in the galactic 
halo, 2006, AIP Conference Proceedings, 861, 807 

[14] J. Diemand, B. Moore, J. Stadcl, Earth-mass dark-matter haloes as the first structures in the 
early Universe, 2005, Nature, 433, 389 

[15] A. M. Green & S. P. Goodwin, On mini-halo encounters with stars, 2007, MNRAS, 375, 1111 

[16] A. Schneider, L. Krauss, B. Moore, Impact of dark matter microhalos on signatures for direct 
and indirect detection, 2010, Phys. Rev. D, 82, 063525 

[17] S. M. Koushiappas, Proper motion of gamma-rays from microhalo sources, 2006, Phys. Rev. 
Lett., 97, 191301 

[18] S. M. Koushiappas, The detection of sub- solar mass dark matter halos, 2009, New J. Phys., 11, 
105012 

[19] M. Kamionkowski, S. M. Koushiappas, M. Kuhlen, Galactic Substructure and Dark Matter 
Annihilation in the Milky Way Halo, 2010, Phys. Rev. D, 81, 043532 



- 11 - 



[20] V. Springel, J. Wang, M. Vogelsberger, A. Ludlow, A. Jenkins, A. Helmi, J. F. Navarro, C. S. 
Frenk, S. D. M. White, The Aquarius Project: the subhaloes of galactic haloes, 2008, MNRAS, 
391, 1685 

[21] J. Stadel, D. Potter, B. Moore, J. Diemand, P. Madau, M. Zemp, M. Kuhlen, V. Quilis, 

Quantifying the heart of darkness with GHALO - a multihillion particle simulation of a galactic 
halo, 2009, MNRAS, 391, L21 

[22] J. F. Navarro, E. Hayashi, C. Power, A. R. Jenkins, C. S. Frenk, S. D. M. White, V. Springel, 
J. Stadcl & T. Quinn, The inner structure of ?CDM haloes - III. Universality and asymptotic 

slopes, 2004, MNRAS, 349, 1039 

[23] J. Diemand, B. Moore, J. Stadel, Convergence and scatter of cluster density profiles, 2004, 
MNRAS, 353, 624 

[24] J. Diemand, M. Kuhlen, P. Madau, M. Zemp, B. Moore, D. Potter & J. Stadel, Clumps and 

streams in the local dark matter distribution, 2008, Nature, 454, 735 

[25] J. F. Navarro, C. S. Frenk, S. D. M. White, , 1996, ApJ, 580, L17 

[26] J. Bullock, T. S. Kolatt, Y. Sigad, R. S. Somerville, A. V. Kravtsov, A. A. Klypin, J. R. 

Primack, A. Dekel, Profiles of dark haloes: evolution, scatter and environment, 2001, MNRAS, 
321, 559 

[27] M. Kuhlen, L. E. Strigari, A. R. Zcntncr, J. S. Bullock, J. R. Primack, Dark energy and dark 

matter haloes, 2005, MNRAS, 375, 387 

[28] A. V. Maccio, A. A. Dutton, F. C. van den Bosch, Concentration, spin and shape of dark 

matter haloes as a function of the cosmological model: WMAPl, WMAP3 and WMAP5 results, 
2008, MNRAS, 391, 1940 

[29] F. Prada, A. A. Klypin, A. J. Cucsta, J. E. Betancort-Rijo, J. Primack, Halo concentrations in 

the standard A cold dark matter cosmology, 2012, MNRAS, 423, 3018 

[30] T. Ishiyama, J. Makino, T. Ebisuzaki, Camma-ray Signal from Earth-mass Dark Matter 
Microhalos, 2010, ApJL, 723, L195 

[31] J. Stadel, Cosmological N-hody simulations and their analysis, 2001, PhD Thesis, University of 
Washington 

[32] E. Bertschinger, Multiscale Gaussian Random Fields and Their Application to Cosmological 
Simulations, 2001, ApJS, 137, 1 

[33] E. Komatsu, Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Cosmological Interpretation, WMAP collaboration, 2011, ApJS, 192, 18 

[34] F. H. Stillinger, Rigorous Basis of the Frenkel-Band Theory of Association Equilibrium, 1963, 
Chem. Phys. 38, 1486 

[35] M. Davis, G. Efstathiou, C. S. Frenk, S. D. M. White, The evolution of large-scale structure in 
a universe dominated by cold dark matter 1985, ApJ, 292, 371 

[36] E. Bertschinger, Self-similar secondary infall and accretion in an Einstein-de Sitter universe, 
1985, ApJS, 58, 39 

[37] J. A. Fillmore, P. Goldreich, Self-similar gravitational collapse in an expanding universe, 1984, 
ApJ, 281, 1 

[38] R. Mohayaee, S. F. Shandarin, Gravitational cooling and density profile near caustics in 
collisionless dark matter haloes, 2006, MNRAS, 366, 1217 

[39] R. Mohayaee, S. F. Shandarin, J. Silk, Dark matter caustics and the enhancement of 

self-annihilation flux, 2007, JCAP, 05, 015 

[40] J. Diemand, M. Kuhlen, Infall Caustics in Dark Matter Halos?, 2008, ApJ, 680, L25 



- 12 - 



[41] M. A Sanchez-Conde, M. Cannoni, F. Zandanel, M. E. Gomez, F. Prada, Dark matter searches 
with Cherenkov telescopes: nearby dwarf galaxies or local galaxy clusters?, JCAP, 12 (2011), 11 

[42] J. Diemand, M. Kuhlen, P. Madau, Dark Matter Substructure and Gamma-Ray Annihilation in 
the Milky Way Halo, 2007, ApJ, 657, 262 

[43] V. Springel, S. D. M. White, C. S. Frenk et al.. Prospects for detecting supersymmetric dark 

matter in the Galactic halo, 2008, Nature, 456, 73 

[44] M. Kuhlen, M. Vogelsberger, R. Angulo, Numerical simulations of the dark universe: State of 
the art and the next decade, 2012, Physics of the Dark Universe, 1, 50 

[45] M. Kuhlen, J. Diemand, P. Madau, The Dark Matter Annihilation Signal from Galactic 

Substructure: Predictions for GLAST, 2008, ApJ, 686, 262 

[46] L. Strigari, S. M. Koushiappas, J. S. Bullock, M. Kaplinghat, Precise constraints on the dark 
matter content of Milky Way dwarf galaxies for gamma-ray experiments, 2007, Phys. Rev. D, 
75, 083526 

[47] L. Pieri, J. Lavalle, G. Bcrtonc, E. Branchini, Implications of High-Resolution Simulations on 
Indirect Dark Matter Searches, 2011, Phys. Rev. D, 023518 

[48] L. Pieri, G. Bertone, E. Branchini, Dark matter annihilation in substructures revised, 2010, 
MNRAS, 384, 1627 

[49] M. A. Sanchez-Conde, F. Prada, On natal concentrations and their implications for dark 
matter annihilation substructure boosts, 2013, in prep. 

[50] F. C. van den Bosch, G. Tormen, C. GiocoH, The mass function and average mass-loss rate of 

dark matter subhaloes, MNRAS, 2005, 359, 1029 

[51] A. Zentner, J. S. Bullock, Halo Substructure and the Power Spectrum, 2003, ApJ, 598, 49 

[52] J. Binney & S. Tremaine, Galactic Dynamics, Princeton University Press (2008) 

[53] C. Giocoli, L. Pieri, G. Tormen, J. Moreno, A merger tree with microsolar mass resolution: 
application to j-ray emission from subhalo population, 2009, MNRAS, 395, 1620 

[54] J. Diemand, M. Kuhlen, P. Madau, Early Supersymmetric Gold Dark Matter Substructure, 
ApJ, 649, 1 

[55] L. Gao, C. S. Frenk, A. Jenkins, V. Springel, S. D. M. White, Where will supersymmetric dark 
matter first be seen?, 2012, MNRAS, 419, 1721 

[56] D. Merritt, A. W. Graham, B. Moore, J. Diemand, B. Terzic, Empirical Models for Dark 
Matter Halos. I. Nonparametric Construction of Density Profiles and Gomparison with 

Parametric Models, ApJ, 132, 2685 

[57] M. Kuhlen, The Dark Matter Annihilation Signal from Dwarf Galaxies and Subhalos, 2010, 
Advances in Astronomy, 162083 

[58] J. Diemand & B. Moore, The Structure and Evolution of Gold Dark Matter Halos, 2011, 

Advanced Science Letters, 4, 297 

[59] A. Klypin, A. V. Kravtsov, J. S. Bullock, J. R. Primack, Resolving the Structure of Gold Dark 
Matter Halos, 2001, ApJ, 554, 903 

[60] M. Ricotti, Dependence of the inner dark matter profile on the halo mass, 2003, MNRAS, 344, 
1237 

[61] A. Tasitsiomi, A. V. Kravtsov, S. Gottlober, A. A. Klypin, Density Profiles of AGDM Clusters, 
2004, ApJ, 607, 125 

[62] J. Diemand, M. Kuhlen, P. Madau, Formation and Evolution of Galaxy Dark Matter Halos and 
Their Substructure, 2007, ApJ, 667, 859 



- 13 - 



[63] W. Dehnen, Phase-space mixing and the merging of cusps, 2005, MNRAS, 360, 892 

[64] M. Zemp, B. Moore, J. Stadel, C. M. Carollo, P. Madau, Multimass spherical structure models 
for N-body simulations, 2008, MNRAS, 386, 1543 

[65] T. Ishiyama, K. Nitadori, J. Makino, 4-45 Pflops Astrophysical N-Body Simulation on K 
computer - The Gravitational Trillion-Body Problem, 2012, eprint arXiv: 121 1.4406 



- 14 - 



